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ABSTRACT 

The epoch of reionization and formation of first stars are inter-linked topics that are of con- 
siderable interest. We use a simplified approach for studying formation of stars in collapsed 
haloes, and the resulting ionization of the inter-galactic medium (IGM). We consider a set of 
ACDM models allowed by observations of CMB temperature and polarization anisotropics 
for this study. We constrain parameters related to star formation with the help of observations. 
We constrain subsets of these parameters independently by using the observed metallicity of 
the inter-galactic medium at z ~ 5 and the requirement that the Thomson scattering optical 
depth due to an ionized IGM as determined for the model from CMB observations be repro- 
duced. We consider a range of initial metalicities for star forming gas, and some variations 
of the initial mass function of stars. We find that a "normal" initial mass function (IMF) may 
satisfy these two constraints with a raised efficiency of star formation as compared to that 
seen in the local universe. Observations require a significant fraction of metals to escape from 
haloes to the IGM. We can also place constraints on the ratio of escape fraction for metals 
and ionizing photons, we find that this ratio is of order unity for most models. This highlights 
the importance of using the constraints arising from enrichment of the inter-galactic medium. 
Ultra-high mass stars or AGNs may not simplify models of reionization in that these may 
produce more ionizing photons but these do not contribute to production of metals and hence 
these help in reducing only the escape fraction for ionizing photons. However, suppression of 
very low mass stars is helpful in that it increases the production of metals as well as ionizing 
photons and hence leads to a reduction in both escape fractions. Such a change is also war- 
ranted by observations of metal poor halo stars in the Galaxy. We also discuss correlations in 
parameters like the efficiency of star formation and the two escape fractions with cosmological 
parameters. 
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1 INTRODUCTION 

Large scale structures in the universe like galaxies and clusters of 
galaxies are believed to have formed by gravitational amplifica- 
tion of small perturbations (Peebles 1980; Shandarin & Zeldovich 
1989; Peacock 1999; Padmanabhan 2002; Bernardeau et al. 2002). 
Much of the matter in galaxies and clusters of galaxies is the so 
called dark matter that is believed to be weakly interacting and 
non-relativistic (Trimble 1987; Komatsu et al. 2008). Dark matter 
responds mainly to gravitational forces, and by virtue of larger den- 
sity than baryonic matter, assembly of matter into haloes and large 
scale structure is driven by gravitational instability of initial per- 
turbations. Galaxies are believed to form when gas in highly over- 
dense haloes cools and collapses to form stars in significant num- 
bers (Hoyle 1953; Rees & Ostriker 1977; Silk 1977; Binney 1977). 
The formation of first stars (McKee & Ostriker 2007; Zinnecker & 



Yorke 2007; Bromm & Larson 2004) in turn leads to emission of 
UV radiation that starts to ionize the inter-galactic medium (IGM). 
The period of transition of the IGM from a completely neutral to 
a completely ionized state is known as the epoch of reionization 
(EoR), e.g., see Loeb & Barkana (2001). 

The study of EoR has been an active area of research in re- 
cent years. Theoretical ideas about the reionization history have 
been constrained by a variety of observations (Fan, Carilli, & Keat- 
ing 2006). For example, observations of Gunn-Peterson troughs in 
AGN spectra at z ~ 6 (Becker et al. 2001; Fan et al. 2006) indi- 
cate that the process of reionization was nearly complete by that 
redshift. Bounds on luminosity function of Lya galaxies at high 
redshifts (Malhotra & Rhoads 2004; Stern et al. 2005; Bouwens et 
al. 2008) also constrain the EoR. These bounds are consistent with 
the conclusion that the IGM was completely ionized by z ~ 6. 
Furthermore, Thomson scattering by free electrons in the IGM of 
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the anisotropic photon distribution that constitutes the CMB leaves 
a signature in its temperature and polarization anisotropy. The op- 
tical depth due to this scattering (Spergel et al. 2003; Dunkley et 
al. 2008) is another probe of EoR. WMAP five-year data indicate 
a value of r ~ 0.084 ± 0.016 that corresponds to z ~ 10 as 
the redshift for instantaneous reionization. Realistic scenarios how- 
ever predict a protracted EoR where the process of ionization starts 
around z ~ 15 and ends by z ~ 6. 

Many possible sources of ionizing radiation have been consid- 
ered, although stellar sources are believed to be the most plausible 
candidates (Yan & Windhorst 2004; Schneider et al. 2006; Choud- 
hury & Ferrara 2007). The AGN density goes down more rapidly 
for 2 > 3 than the density of star-forming galaxies (Miralda- 
Escude & Ostriker 1990; Madau, Haardt, & Rees 1999; Haehnelt et 
al. 2001) and therefore AGNs are not expected to contribute signif- 
icantly to the ionizing radiation. X-rays from low mass quasars, X- 
ray binaries and supernova remnants are constrained by the soft x- 
ray background observed today (Dijkstra, Haiman, & Loeb 2004). 
Particle decays can also play only a minor role in reionization 
(Bharadwaj & Sethi 1998; Mapelli & Ferrara 2005). In the present 
work, we assume that it was radiation from early stars that ionized 
the IGM and ignore other possibilities. 

The total photon emissivity of early stars is poorly known. 
Studies of reionization typically use observations with semi- 
analytic models of heating and ionization of the IGM where ef- 
ficiency of star formation, evolution of star formation rates, the 
number of ionizing photons emitted per baryon in stars, etc. are 
parameterized in some manner (Chiu & Ostriker 2000). Given the 
complexity of most of these approaches, and the number of param- 
eters, it is often impractical to scan the parameter space. The main 
lesson we learn from these studies is that we may not require ex- 
traordinary physical processes in order to satisfy available observa- 
tional constraints of reionization. We make an attempt to simplify 
modeling of star formation and other astrophysical aspects of the 
problem. This allows us to reduce the number of free parameters 
in this sector while retaining many significant astrophysical rela- 
tionships. Statements can then be made about quantities like the 
escape fraction of UV photons, /osc, 7 and their correlations with 
the cosmological parameters. 

Large uncertainties exist in parameters related to early star 
formation, e.g. the escape fraction of UV photons, f asc ,-y, the ini- 
tial mass function (IMF) and the efficiency of star formation, /* 
(Bunker et al. 2004). It has been suggested in the literature that a 
top-heavy IMF with very massive stars is not necessarily favored to 
satisfy the reionization and metal enrichment constraints (Daigne et 
al. 2006; Venkatesan & Truran 2003). 

In this work we assume that early star formation happens pre- 
dominantly during formation and major mergers of haloes, and oc- 
curs as a short lived burst. Photon emissivities and metal yields can 
then be calculated using population synthesis models for different 
IMFs (Leitherer et al. 1999; Bruzual & Chariot 2003). We then test 
if these scenarios generate enough photons to ionize the universe 
by z ~ 6 by comparing with the observed Thomson scattering op- 
tical depth. We also require the models to satisfy constraints arising 
from the observations of the metal content of the IGM. This con- 
strains the amount of processed elements that escape from the ISM 
to the IGM. We can use models for outflows as a guide and put 
constraints on the efficiency of star formation, or use "reasonable" 
values of the efficiency of star formation to constrain the metal es- 
cape fraction. We can also combine the two constraints to scale out 
efficiency of star formation. The questions we address in this work 
are: 



• To combine constraints of enrichment of the inter-galactic 
medium with observations of reionization, and check whether the 
extra information can provide constraints on the initial mass func- 
tion during early star formation. 

• Can the combined constraints be used to make useful state- 
ments with regard to other potential sources of ionization? 

• Is there any correlation between the parameters that describe 
star formation and cosmological parameters? 

Relevant observations are described in §2. The reionization 
model that we use is discussed in §3. Results and discussion appear 
in 54 and 55. 



2 OBSERVATIONS AND MODELS 

The IGM occupies most of the space and a substantial amount of 
matter in the Universe. Indeed, it is believed that at least half of the 
baryons in the universe are in the IGM. Thus it is not surprising that 
the observations of IGM dominate when we discuss constraints on 
models of reionization. We introduce observations that are used for 
constraining models in this paper. 

Observational constraints that we use are essentially two: 
metallicity of the IGM as seen in quasar absorption systems, and 
the Thomson scattering optical depth from the cosmic microwave 
background. A third observation that can potentially be used is that 
of the cosmic stellar matter density, which we model in a man- 
ner explained below. We do not discuss constraints arising from 
this type of observations as the observational bounds on models 
are not very strong at present. It is not possible to use other ob- 
servational constraints (Choudhury & Ferrara 2007) in the global 
averaged model discussed here. 

Observations of galaxies at high redshifts can be used to infer 
the density of stellar matter, p», at those redshifts. One way of esti- 
mating this quantity is to use the luminosity function of galaxies in 
various wavelength passbands and combine these with population 
synthesis models and an assumed IMF (Madau et al. 1996; Lilly 
et al. 1996; Bouwens et al. 2007). We do not use these observa- 
tions here as current observations do not provide sufficiently strong 
constraints. This is expected to change in coming years with better 
observations. 

Instead in our work we assume that stars form only in viri- 
alised haloes. We take the minimum mass of star-forming haloes 
to be 10 s M©, this is a proxy for haloes with a virial temperature 
of 10 4 K. We do not take into account star formation in haloes 
of around 10 6 Mq that is aided by molecular cooling (Tegmark 
et al. 1997). We also ignore the effects of feedback that raise the 
Jeans mass to around 10 9 M '. We can then obtain metal and pho- 
ton yields by using population synthesis models: we use the STAR- 
BURST99 code (Leitherer et al. 1999; Vazquez & Leitherer 2005). 
This assumption allows us to connect the rate of change of stellar 
mass to the rate of change of the total mass contained in massive 
haloes. We obtain this rate from the Press-Schechter formalism for 
a Gaussian PDF as 



F(m,z) = 



2 i<5 c (l + . 



it a(m)d + (z) ^ \ v 2u 2 (m) J <91oga 
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-S 



d log d+ 



(1) 

where an overdot denotes derivative with respect to the cosmic 
time, a prime denotes derivative with respect to the redshift, and 

1 We have checked that including the effects of radiative feedback in- 
creases the required star formation efficiency by around 20%. 
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0.0004 
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Kroupa 


0.5 


100 


0.0004 


9280 


0.0167 


50.0 


3. 


Kroupa 


0.1 


100 


0.001 


6297 


0.0159 


50.0 


4. 


Salpeter 


1 


100 


0.001 


11237 


0.0283 


50.0 


5. 


Kroupa 


0.1 


100 


0.02 


3996 


0.0261 


50.2 



Table 1. Various IMFs used in this study are summarized here. Mi ow and M b j gb are the lower and upper mass cut-offs for the IMFs respectively. Z- lnput 
denotes the metallicity of the gas from which stars are formed. Nj is the number of ionizing photons produced per baryon in stars and p is the metal yield per 
baryon in them. 7V 7 and p are obtained using population synthesis models. The last column lists the value of Af 7 / eS c,7/* for the WMAP5 best-fit model with 
the corresponding IMF. This quantity is proportional to the number of photons escaping into the IGM for every baryon inside a collapse halo. 



F is the fraction of haloes with mass greater than m (Press & 
Schechter 1974). The critical density for spherical collapse is sym- 
bolized by S c , and a 2 (m) is the variance in the initial density fluc- 
tuation field when smoothed with a top-hat filter of a scale corre- 
sponding to mass m. The rate of growth for perturbations in the 
linear theory is denoted by d + (z). 

The total amount of baryons added to haloes of mass greater 
than 10 8 M Q is taken to be the amount of gas available for star for- 
mation. Gas already present in haloes is not considered for star for- 
mation. We do not consider the contribution of minihaloes as these 
do not contribute significantly to the total star formation due to ra- 
diative feedback (Trenti & Stiavelli 2009). We define the efficiency 
of star formation, /» as the fraction of this gas that is converted into 
stars. The rate of change of stellar mass can now be written as 

p«(z) = fM hP cF{W 8 M Q ,z), (2) 

where p c is the critical density and fib is the density parameter for 
baryons. Here we have ignored mass lost by stars through winds, 
outflows, and supernovae. This can be taken into account using, 
for example, population synthesis models 2 . Equation (2) illustrates 
the small number of parameters and approximations that go into 
estimating p\ in our formulation. 

Observations of absorption systems in quasar spectra have 
been used to put constraints on the average metallicity of the IGM 
(Cowie et al. 1995; Songaila 1997; Ellison et al. 2000; Simcoe, Sar- 
gent, & Rauch 2004; Becker, Rauch, & Sargent 2008; Ryan- Weber 
et al. 2009). These observations indicate that the amount of C IV 
in the IGM does not evolve significantly between 2 < z < 5.5 
(Songaila 2001; Becker, Rauch, & Sargent 2008; Ryan- Weber et al. 
2009). It is not yet clear whether the IGM has been contaminated 
by metals throughout, or if the enriched regions of the IGM are re- 
stricted to the neighborhood of galaxies and filaments. There are 
also issues related with understanding the ionization state to map 
from absorption by a particular species to metallicity (Schaye et 
al. 2003). Observations also support a correlation between density 
and metallicity of the IGM, indicating that regions in proximity of 
galaxies are enriched to higher level than regions of IGM far away 
from galaxies (Schaye et al. 2003; Pieri, Schaye, & Aguirre 2006; 
Scannapieco et al. 2006). We assume that the IGM is uniformly 
enriched at the level indicated by Songaila (2001) at z ~ 5.5. 

Winds and outflows are expected to be the dominant processes 
that lead to ejection of metals from the inter-stellar medium (ISM) 
of galaxies. There is considerable evidence in favor of this mech- 
anism as observations have detected outflows around almost all 

2 A starburst with a Kroupa IMF and an initial metallicity of 0.02 loses 
about 10% of its mass to the ISM through these effects under normal as- 
sumptions. 



galaxies at high redshifts (Pettini et al. 2001; Frye, Broadhurst, & 
Bem'tez 2002). The fraction of processed metals that can be de- 
posited from the ISM to the IGM without disturbing the IGM in an 
observable manner is not known. Several authors often assume that 
around 1% of metals produced in galaxies can be ejected and de- 
posited in the IGM. This may also be computed from first principles 
in detailed models (Daigne et al. 2004, 2006; Samui, Subramanian, 
& Srianand 2008). In these models supernova-driven outflows are 
responsible for the IGM enrichment. The efficiency of these out- 
flows depends on the star formation efficiency, the IMF and the ef- 
ficiency of winds (Madau, Ferrara, & Rees 2001; Scannapieco, Fer- 
rara, & Madau 2002; Scannapieco 2005; Furlanetto & Loeb 2003). 
For a star formation efficiency of 10% the volume filling factor of 
the ejecta can be 20-30% and at z ~ 3 the IGM metallicity could 
be around [-3] as detected by Songaila (2001). 

The amount of metals produced per baryon in stars can be 
computed once we fix the initial mass function (IMF) of stars, and 
metallicity of the star-forming gas (Leitherer et al. 1999; Vazquez 
& Leitherer 2005). We can write 

nz = /„/c SC ,zfib — F(10 8 Mq,z)p (3) 
Trip 

for the number density of metals that reaches the IGM. Here / csc ,z 
is the fraction of total metals produced that is deposited in the IGM 
and p is the metal yield of the stars per baryon. Note that we as- 
sume the same escape fraction for metals from galaxies of different 
masses, whereas it is far more likely that low mass galaxies lose 
nearly all the metals from the ISM and more massive galaxies lose 
very little (Dekel & Silk 1986). 

Finally, observations of the temperature and polarization 
anisotropies in the CMB provide a constraint on the EoR. Free 
electrons produced during reionization scatter the CMB photons 
and suppress temperature anisotropies on scales smaller than the 
Hubble radius at that time. The suppression is of the form Cf = 
e- 2T C[ where 

r = j n e (t) a T cdt (4) 

is the Thomson scattering optical depth. Here n e (t) is the number 
density of free electrons and ot is the Thomson scattering cross- 
section. This damping is degenerate with the amplitude of the pri- 
mordial power spectrum. The degeneracy is broken by detection of 
a polarization anisotropy for scales greater than the Hubble radius 
at the reionization redshift, an effect that dominates at scale of the 
Hubble radius at EoR and has an amplitude proportional to r. Any 
model of reionization must reproduce the observed value of optical 
depth. For a forecast on anticipated improvements of observational 
estimation of r, please see (Colombo & Pierpaoli 2009). 
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3 REIONIZATION MODEL 

The reionization history depends on the star formation history of 
the universe, which in the simplified models is closely related to 
the halo formation history. The IMF of stars and the escape frac- 
tion for ionizing photons then give us the number of ionizing pho- 
tons that are available as a function of time. These can then be used 
to compute the evolution of the neutral or ionized fraction of gas 
in the universe. As mentioned above, we assume that star forma- 
tion is triggered during formation of haloes. As most time scales 
of interest are longer than the dynamical time scale over which the 
bulk of star formation takes place, we assume star formation to be 
instantaneous in our model. 

We consider a global averaged evolution of ionized fraction 
instead of following evolution of HII regions around haloes, the 
approach used in most studies (Chiu & Ostriker 2000; Sethi 2005). 
Further, we assume that during reionization, a region is either neu- 
tral or completely ionized. With these assumptions, the evolution 
of the ionized fraction evolves as: 

x = — aBC 2 nux + <j p ynu_c{l — x) (5) 
y = -a p ynHc(l-x) + m p f*f eS c, 7 FN 7 . (6) 

Here x is the fractional volume that is ionized, and y is the 
number of ionizing photons per baryon. a p denotes the effective 
cross-section of photoionization, as is the recombination coeffi- 
cient for all levels except the ground state of neutral hydrogen, and 
m p denotes the mass of a proton. 

The first term on the right hand side of equation (5) de- 
scribes recombination. C is the clumping factor defined as C = 
{n 2 H ) / {uh) 2 ■ This term usually involves square of the ionized frac- 
tion but in our model we assume that the ionized fraction is either 
unity or zero. This, when used in volume averaging over the uni- 
verse with an additional assumption that the clumping is the same 
in ionized and neutral regions, leads to a linear dependence. In the 
process of averaging, the meaning of x changes from the ionized 
fraction to the volume filling fraction of the ionized regions. We 
can express this in terms of equations: 

1 / nWdV =<!p.J x 2 dV =^lj xdV = {n 2 H )x. (7) 

We have assumed that the clumping factor is the same in all parts of 
the universe, this allows us to take (n H ) outside the integral. The 
third equality in equation (7) follows from the definition of x as a 
filling fraction. We also expect C to change with the evolution of 
clustering. We take this dependence to be of the form (Iliev et al. 
2007) 

C 2 ~ 26.2917 cxp [-0.18222 + 0.003505z 2 ] . (8) 

Sources of ionizing radiation are represented in the last term 
of equation (6), F being related to the formation rate of collapsed 
haloes. This is obtained from the Press-Schechter formalism as de- 
scribed above. denotes the number of photons produced per unit 
mass of star formation. Ionization of neutral hydrogen is described 
in the last term on the right hand side of equation (5). This term 
occurs in both equations. We neglect the contribution of collisional 
ionization. 

We solve these equations numerically for different cosmolog- 
ical models. The system of equations (5) and (6) is "stiff," since 
(x — 0, y = 0) is a stable point and time scales for evolution of 
x and y are very different. Further, x is bounded from above (by 
unity) while y is not. Thus the usual forward differencing methods 
do not give accurate solutions easily. We bypass this problem by 



noting that during the process of reionization almost every ionizing 
photon will be immediately absorbed by the medium 3 . This means 
that the two terms in the right hand side of the second equation are 
of the same order till x becomes nearly equal to 1, whereas the left 
hand side is much smaller and may be assumed to be zero. This re- 
duces the system of equations to a single equation, which can now 
be solved using forward differencing methods. Note that this ap- 
proximation is not valid when x approaches 1, although in practice 
the approximate solution is fairly accurate up to x ~ 0.9. Indeed, 
if we use the approximation up to a; = 1.0 then we make an error 
in estimation of r of less than 5%. 

We take aa = 1.0 x 10~ 13 cm 3 sec" 1 , ignoring its depen- 
dence on temperature. This dependence is fairly weak at tempera- 
tures of interest. We use a p — 6.30 x 10 _ cm . We thus assume 
that most of the ionizing radiation is around the Lyman limit. The 
number of ionizing photons released per baryon of stars formed, 
denoted by , depends on the initial mass function (IMF) of the 
stars. We obtain this number from the STARBURST99 stellar popu- 
lation synthesis code 4 (Leitherer et al. 1999; Vazquez & Leitherer 
2005). 



4 DISCUSSION 

Our aim here is to study a variety of models with varying cos- 
mological parameters as well as parameters related to star forma- 
tion and enrichment. We consider a random subset of flat ACDM 
models allowed by WMAP5 (Komatsu et al. 2008; Dunkley et 
al. 2008). We do not consider models with massive neutrinos or 
a non-vanishing tensor component, or models where the primor- 
dial power spectrum deviates from a pure power law. We use only 
WMAP constraints for limiting cosmological parameters. We used 
the MCMC chains made available by the WMAP team (http: 
/ /lambda . gsf c . nasa . gov/). We considered a random sub- 
set of all models allowed with a confidence level of 68% from the 
MCMC chains. We studied a handful of models for parameters re- 
lated to star formation; details of these are given in Table (1). The 
table lists the IMFs used in our study 5 . We have also listed the 
amount of ionizing photons produced per baryon in stars, and the 
total metal yield per baryon for these IMFs. These numbers also 
depend on the metallicity of gas from which stars form and this is 
listed in the table as well. It is interesting to note that for a given 
IMF, as the metallicity increases, the production of ionizing pho- 
tons per baryon comes down but the amount of enriched material 
returned to the ISM increases. Gas that forms the first stars is likely 
to have primordial abundance (Olive, Steigman, & Walker 2000). 
In our analysis of all the models, we keep metallicity fixed and 
hence it is appropriate to use low values of input metallicity. 

3 There are two approximations being discussed here: one approximation 
is that all UV photons are available for ionizing atoms. This is not really 
true as ionized regions can be expected to host an ionizing background. The 
other approximation is in solving the equations that we arrive at with the 
first approximation. 

4 We consider instantaneous starbursts with a fixed total stelar mass of 10 6 
Mq. We use the Geneva evolutionary tracks for models with metallicity 
0.001 and the Padova tracks with AGB stars for models with metallicity 
0.0004. 

5 We should note that there is considerable uncertainty in the shape of 
the IMF in the local neighborhood (Kroupa 2002; Conroy, Gunn, & White 
2008). This can easily have a significant impact on our conclusions. The 
uncertainties introduced by other assumptions and approximations should 
be seen with the uncertainty in the IMF as the reference. 



Constraints on Early Star Formation 5 




0.96 0.98 



Figure 1. The products /„ / esc ,z (top row) and /* / eS c,7 (bottom row) against cosmological parameters for model 1 of Table (1). The star symbol denotes the 
WMAP5 best fit model. These points also represent lower bounds on / esc ,z and /esc, 7. See text for details. 



Metallicity of the IGM constrains the product /*/ eS c,z for a 
given model. Similarly we constrain the product /*/ C sc,7 with the 
optical depth due to reionization for the CMB. Figure (1) shows 
these products for all cosmological models studies here, when the 
star formation parameters for model 1 in Table (1) are used. Given 
that the efficiency of star formation can at best be 100%, i.e., /» < 
1, the points in Figure (1) also represent lower bounds on / eS c,z and 
/esc, 7- These are shown as a function of the slope of the primordial 
power spectrum (n s ), optical depth due to reionization (t C mb), and, 
amplitude of clustering at the scale of 8 h~ Mpc {as). The best fit 
WMAP5 model is marked in each panel as a star. We find that there 
is some correlation between the lower bound on /» /esc, 7 and t C mb , 
and also between /* / CS c,z and n s . There are weak correlations with 
other cosmological parameters but nothing as remarkable as the 
two mentioned above. 

The efficiency of star formation is likely to be much less than 
unity in any realistic scenario. Indeed, if we try to keep the different 
efficiencies and escape fractions at the same order then we require 
these to be around 0.1 — 0.15, or 10 — 15%. While the escape 
fraction for ionizing photons and star formation efficiency we have 
obtained are comparable to those found in other studies, these are 
higher than the values seen in local galaxies. In particular, it is not 
clear if it is possible to expel 10% of the metals from the ISM to 
distant parts of the IGM using known physical mechanisms. Before 
commenting on the numbers, let us consider the sensitivity of the 
result to our assumptions. 

• If the IMF has a low mass cutoff that is higher than the 0. 1 M 



used for model 1 from Table (1), then a larger fraction of mass goes 
into high mass stars that produce the ionizing photon flux and the 
enriched material. This can lower the required / CS c,z by a signif- 
icant amount. The product /,/ eS c,z for model 2 from Table (1) is 
shown in Figure (2). There is some evidence that there are more 
intermediate mass stars in the population of metal poor stars in the 
halo of the Galaxy as compared to metal rich stars, if we normal- 
ize the two distributions at low stellar masses (Tumlinson 2007; 
Komiya et al. 2007). Thus a higher cutoff for Mi ow may be re- 
quired for explaining other observations. 

• Our analysis assumes that the metallicity of gas that forms 
stars is fixed. If we do a self-consistent analysis where this is al- 
lowed to evolve, the gas metallicity gradually increases. It is then 
clear from Table (1) that later generations of stars will enrich the 
ISM faster. As an illustration, we can see the results of analysis 
with higher fixed input metallicity for models 2-4 in Figure (2). 
Our estimates show that / eS c,z required to satisfy observations can 
come down by a few tens of percents due to this. There is a cor- 
responding increase in / CS c,7 due to stars with higher metallicity 
producing fewer ionizing photons. 

• We have assumed that haloes with mass above 10 s M can 
form stars. The standard approach is to assume that radiative feed- 
back during the EoR increase this by about an order of magnitude in 
regions that have been photo-ionized (Efstathiou 1992). We do not 
take this into account as it has been pointed out that the actual ef- 
fect may only be to reduce the efficiency of star formation in lower 
mass haloes (Mesinger & Dijkstra 2008). If we do consider the ef- 
fect of radiative feedback as disabling star formation then it leads 
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Figure 2. The products /*/ eS c,z and /*/ CS c,7 for model 2 (left column), model 3 (center column) and model 4 (right column) of Table (1). Corresponding 
plots from Figure (1) are superimposed in grey. Filled stars denote values for the WMAP5 best-fit models. These points are represent lower bounds on / esCi z 



and / e , 



v. See text for details. 



to a reduction in total gas available for star formation, and hence 
requires slightly higher efficiencies and escape fractions. We find 
that this effect requires an increase in the two products by about 
20%. 

• We have assumed that the universe is enriched uniformly. It 
may very well happen that the enrichment process is effective only 
in the vicinity of galaxies. In such a case only overdense regions are 
enriched. The escape fraction of metals required can be lowered by 
as much as a factor of two if this is the case. 

Thus we may require only around 5% of the ISM to be ejected to 
the IGM on an average in models with /» ~ 0.2 This is compa- 
rable with semi-analytic ab initio models of early star formation, 
outflows and IGM enrichment that have been studied in the litera- 
ture. 

We have assumed the same loss fraction for ISM for galax- 
ies over the entire range of masses. This, of course, is not true. We 
expect that the low mass galaxies can potentially disperse a large 
fraction of the ISM in supernova explosions but larger galaxies can 
retain most of their ISM (Larson 1974; Dekel & Silk 1986). If most 
of the IGM enrichment is done by metals that form in dwarf galax- 
ies then the constraint is not very stringent. In most models, the 
fraction of mass in galaxies with a halo mass of less than 10 10 Mq 
is larger than 10% even at z ~ 6. If these galaxies lose a significant 
fraction of the ISM on an average and heavier galaxies lose very 
little mass then we can comfortably satisfy the constraints from en- 
richment of the IGM. 

We now turn our attention to the ratio of escape fractions for 




Figure 3. Ratio / OS c,z//csc,7 for the first IMF listed in Table (1). See text 
for details. 



photons and metals. Figure (3) shows the ratio / CS c,z//csc,7 for the 
first IMF listed in Table (1). It is interesting to note that the ratio 
/csc,z//csc,7 is of order unity, differing from unity by at most a 
factor of a few. Thus the fraction of ionizing photons that escape 
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galaxies is broadly of the same order as the fraction of metals that 
must leak into the IGM in order to explain the observed enrichment 
of the IGM. We have not plotted this ratio for other IMFs in Ta- 
ble (1) as the expected change can be seen from Figure 2. Indeed, 
the change in the ratio is less than a factor two as we consider IMFs 
listed in rows 2 — 4 of the Table (1). The last row in the Table (1) 
is more appropriate for late time star formation and we need not 
discuss that here. 

Most studies of reionization have tended to focus on the es- 
cape of ionizing photons. In order to satisfy observations of r or the 
luminosity function of high redshift galaxies, these often invoke a 
top heavy IMF for the first generation of stars (Cen 2003; Haiman 
& Holder 2003; Wyithe & Loeb 2003; Bromm 2004). Other sources 
of ionizing radiation like AGNs, magnetic fields and decaying dark 
matter particles have also been studied (Pierpaoli 2004; Schleicher, 
Banerjee, & Klessen 2008). It is interesting to note that all such 
modifications lead to an enhanced production of ionizing photons 
without affecting the production of metals in a significant manner. 
This is because the very massive stars are expected to implode and 
do not enrich the ISM with the products of nuclear fusion that takes 
place in the core. The ratio / OS c,z//osc,7 changes on addition of 
extra sources of ionizing radiation. In view of the arguments pre- 
sented above, all modifications that have been discussed so far lead 
to a smaller / eS c,7- In other words, the ratio plotted in Figure (3) 
should be thought of as a lower bound. 

An important implication of this is that the constraints from 
enrichment of the IGM require certain amount of star formation, 
and this requirement needs to be satisfied even when we invoke 
other sources of ionizing radiation during the epoch of reionization. 
That is, adding new potential sources of ionizing radiation can be 
helpful only in lowering the escape fraction of ionizing photons and 
not in lowering the amount of star formation 6 . We may even end up 
with a scenario where a much larger fraction of processed elements 
need to be transferred from the ISM to the IGM as compared to the 
fraction of UV photons escaping from galaxies. 



5 CONCLUSIONS 

We have compared simple models for star formation in the early 
universe with two observational constraints. The simplicity of the 
model allows us to consider variation in cosmological parameters 
as well. We present a summary of our results here: 

• The product of star formation efficiency and escape fraction 
of ionizing photons, /*/ es c,7, is correlated with the optical depth 
due to reionization. 

• The product of star formation efficiency and escape fraction 
of ISM, /*/esc,z, is anti-correlated with the index of the primordial 
power spectrum. 

• These are weak correlations, in the sense that the values for 
the two products do not change strongly for small changes in the 
cosmological parameter in question. 

• We do not find any other correlation amongst parameters of 
star formation and cosmological parameters. 

• We are able to satisfy observational constraints with the stan- 
dard initial mass function for stars observed in the local universe 
(Kroupa 2002), and with reasonable values for star formation ef- 
ficiency and escape fractions for photons and ISM. Given that the 

6 An exception is the scenario where the universe is not enriched through- 
out. In such a case even / ea c,z can be reduced by a significant amount. 



local IMF itself is somewhat ill constrained, this implies that we do 
not require a significant evolution of the IMF in order to explain 
observations considered here. 

• Small variations in the IMF, indicated by observations of metal 
poor stars in the Galaxy, reduce the efficiency of star formation and 
the escape fractions required for the standard IMF. 

• Approximations used by us in the model do not change the 
overall numbers by more than 10 — 20%. Indeed, different approx- 
imations change numbers in different directions so we can consider 
the overall results to be fairly robust. 

• Our model allows us to estimate the ratio of the two escape 
fractions. We find that the two escape fractions are of the same 
order. 

• If we consider other potential sources of ionizing photons then 
the required escape fraction for photons can come down, however 
the escape fraction for processed elements does not change. Indi- 
rectly, the required amount of star formation is required to remain 
the same unless there is some very efficient mechanism for trans- 
porting processed elements into the IGM while keeping the escape 
fraction of ionizing photons low. 

The most important conclusion of this work is that star forma- 
tion without a significant evolution of the IMF is sufficient for sat- 
isfying the two constraints considered here. The escape fractions, 
and/or the star formation efficiency is required to be higher than we 
see in local galaxies. One can consider other sources of ionizing 
radiation, indeed at least some of these must be present. But as we 
have pointed out, these help in reducing only the escape fraction 
for ionizing radiation as none of the other potential sources help 
in transporting enriched material from the inter-stellar medium to 
the inter-galactic medium. This highlights the significance of the 
constraint arising from enrichment of the IGM for epoch of reion- 
ization studies. 
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